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Abstract 

The EDDA-Detector at the Cooler-Synchrotron COSY/Jiihch has been operated 
with an internal CH2 fiber target to measure proton-proton elastic scattering dif- 
ferential cross sections. For the data analysis knowledge of beam parameters, like 
position, width and angle, are indispensable. We have developed a method to obtain 
these values with high precision from the azimuthal and polar angles of the ejectiles 
only, by exploiting the coplanarity of the two final state protons with the beam 
and the kinematic correlation. The formalism is described and results for beam 
parameters obtained during beam acceleration are given. 
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1 Introduction 



In nuclear and particle physics experiment knowledge of the distribution of 
scattering vertices is a crucial ingredient of data analysis. In most cases ex- 
periments are designed to yield enough position information to trace detected 
particles back to the origin. However, even with a much simpler detector, 
events with particularly simple topologies like elastic scattering may allow to 
extract the vertex distribution with much lesser experimental information. 

The EDDA experiment [1,2] at the COSY accelerator [3] in Jiilich is dedicated 
to the measurement of pp elastic scattering over a wide energy (0.5-2.5 GeV) 
and angular range (30° to 90° in the cm.). The experimental program in- 
cludes measurements of the unpolarized differential cross-section [1,2], analyz- 
ing power [4] and spin-correlation parameter [5] . Data taking for unpolarized 
cross section was done only with the outer part of the EDDA-detector, shown 
in Fig. 1, detecting in coincidence a single point of incidence for each of the two 
protons, originating from elastic scattering of the internal proton beam from 
5x4 //m^ thin CH2-fiber targets. However, for the reconstruction of scattering 
angles, detector calibration, and luminosity determination detailed knowledge 
of the beam-target overlap is important. Since measurements were taken dur- 
ing acceleration of the COSY-beam, such that beam position and width are a 
function of beam-momentum, it was mandatory, that these quantities can be 
determined from the scattering data itself. 

Two-body kinematics require the two protons to be coplanar with the beam. 
Thus deviations from coplanarity of proton-proton scattering events with the 
detector symmetry axis bears information on the position and angle of the 
COSY beam axis with respect to the detector. Along these hues we developed 
a method to extract information on beam and target parameters like width, 
position and angle from a single-layered detector described in this paper. In 
sections 2 and 3 we will give a brief description of the experiment and the 
selection of scattering data. Section 4 will explain the formalism and in section 
5 some examples for the application of this method will be discussed. 



2 Experiment 

The EDDA Detector (cf. Fig. 1) measures the point of interception of the two 
final-state particles of proton-proton scattering events with a cylindrical sym- 
metric scintillator array around the COSY beam pipe. It consists of a layer 
of 32 overlapping scintillator bars (B) of triangular cross section detecting 
the azimuthal angle (p and a layer of 2x29 semirings, made either from bulk 
scintillator (R) or scintillating fibers (F), detecting the position z of incidence 
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Fig. 1. Schematic drawing of the EDDA detector 

along the symmetry axis of the detector and bears information on the polar 
angle 6'iab- The angular resolution is improved by about a factor of five with re- 
spect to the granularity by comparing off-line the light output of two adjacent 
detector elements [6]. 

Thin CH2- or carbon-fibers, stretched horizontally on a 30 mm wide fork, are 
used as targets. They can be moved vertically in and out of the beam by 
a magnetic drive. Beam lifetimes of 3 to 10 s, depending on beam energy, 
have been achieved. Data is taken during and/or after acceleration of the 
COSY-beam and accumulated over many machine cycles to obtain statistical 
precision. Hits in opposing bars serve as a coplanarity trigger. In addition, 
elastic scattering kinematics 

cot ^lab,! • cot eiab,2 = ^^2^ = 7c.m.> (l) 

with 7c.ni. being the Lorentz-factor of the cm. in the laboratory and R^i the 
radius of the ring-layer, is tested by a coincidence of two semiring with appro- 
priate locations in z. 



3 Event Selection 

Coplanarity and the kinematic correlation (cf. Eq. (1)) can be tested with 

refined precision offline. Then two prongs detected in coincidence are trans- 
formed to the cm. system, with the hypothesis of being elastically scattered 
protons. For pp elastic scattering events these prongs should be back-to-back 
in the cm. system, and the angular deviation a from this signature, the so- 
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Fig. 2. Distribution of the kinematic deficit a for scattering events at a beam mo- 
mentum of 1.455 GeV/c. The upper hmit for accepted elastic scattering events is 
shown as the vertical line. 

called kinematic deficit, should be small. A typical distribution is shown in 
Fig. 2. The apparent contribution from the carbon contents of the target is 
measured separately using pure C-targets and subtracted statistically. To this 
end two data sets taken with a CH2- and a C-target under the same exper- 
imental conditions are subject to an identical offline analysis. The relative 
normalization of the two samples with respect to proton-carbon reactions is 
derived from the counts in the range 10° to 15°in a (cf. Fig. 2), mainly popu- 
lated by quasi-free scattering on protons bound in carbon. This normalization 
factor is used to subtract the contribution from pC scattering in all acquired 
experimental spectra, so that a clean pp elastic scattering sample is available. 
Monte-Carlo simulations showed this procedure also takes care of the small 
contribution from inelastic pp-reactions in first order. Details of the event 
reconstruction are described elsewhere [1,2]. 



4 Formalism 



In this section we will use a coordinate system (x, y, z) attached to the EDDA 
detector as shown in Fig. 1, with the z-nxis being the symmetry axis of the 
detector pointing downstream from the target (at z = 0) with respect to the 
COSY-beam, the y-axis pointing upward in the lab, and the x-axis pointing 
outboard of the COSY-ring. If not stated otherwise all quantities will be given 
in this detector coordinate system. Widths given as standard deviation are 
denoted by cr, while S signifies FWHM (full width at half maximum) 

The distribution of scattering vertices must be along the fiber target, which 
have negligible (4-5 /xm) spatial extension in the y- and ^-directions. The 
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Fig. 3. Geometry for a beam parallel to the detector symmetry axis z as explained 
in the text. 

longitudinal target position along z can be determined from scattering data 
using Eq. (1) and defines the zero of the coordinate system in this direction. 
Finally, the parameters yet to be determined are 

fiber target position in y at z=0. 
Xft beam position in x at 2;=0 and y^. 

x'^ horizontal beam angle (dx/ds) at z—^. 

yl vertical beam angle (dy/ds) at z—0. 

Sxb^FWHM size of the beam in x at z—0 and yb. 
5(f) detector resolution in the azimuthal angles and 02- 

For simplicity the formalism will first be described for a beam parallel to the 
detector symmetry axis, i.e. a;|,= y[— 0. The generalization to tilted beams 
will be the topic of Section 4.2 

J,..! Beam Parallel to the Detector Axis 

Consider a beam displaced by {xb,yb) in the x,y-plane from the detector sym- 
metry axis. In a projection on the x,y-plane of an elastic scattering event 
(prongs pi and p2 in Fig. 3) a beam displacement leads to a deviation from 
coplanarity when the vertex is assumed to be at {x,y)—{0,0) (dotted lines). 

The EDDA detector measures the point of interception of these two prongs 
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with the bar layer at a mean radius of 7? = 164 mm. Connecting these two 
points by a straight hne its distance d to the origin of the coordinate system 
is a function of the beam position. If we now define an axis d perpendicular 

— * 

to this hne (fat dashed hne in Fig. 3) we can calculate the angle 4>d of the d 
axis to the x-axis 

(Pd = — 7, (2) 



and the position d along this axis: 

d — —Rcos — - — (3) 



— * 

The orientation of the d axis has been chosen to point to the "left" , to the effect 
that (t)d e [— 7r/2,7r/2] must be forced by adding or subtracting tt from the 
expression given in Eq. (2) (For the example of Fig. 3 vr must be subtracted). 
The negative sign in Eq. (3) is due to this choice. Eqs. (2) and (3) are an 
eventwise transformation from the observables {4>i,4>2) to {d,4>d)- 

The connection between d and the beam position is evident. Using basic 
trigonometric relations (for a derivation see Appendix A) - or by a simple 
geometric derivation using Fig. 3 - one obtains: 

d^Xb cos 0d + yb sin (pa (4) 



Now, the idea is to consider the d-distribution of events with different 0^ in 
order to deduce Xb and y^. 



4.1.1 Experimental Input and Analysis 

For each event (pd and d, are calculated. The events are sorted into spectra 
vs. d for different ranges in (f)d. Events with \(f)d\ < 10° are discarded, since 
here the readout of the semirings causes trigger inefficiencies. A (i-distribution 
is accumulated for 16 ^^-intervals from -90° to -10° and 10° to 90°, each 5° 
wide. These 16 spectra, an example is shown in Fig. 4 (a), are the input to the 
analysis. For investigation of the time dependence of beam parameters this 
set of spectra is accumulated for small time intervals, usually 25, 50 or 100 ms 
wide, separately. 

The experimental ^-distributions can be fitted by a Gaussian 

<d) = r— exp -2 — (5) 

V27r(Td V J 
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Fig. 4. (a) Variation of the d-distribution for different bins in 4>d- The fit (sohd 
Gaussian) to the data (crosses) yields a beam displacement of Xb = —2.51 mm, 
Hf, = —0.41 mm. The mean value is shown as the solid vertical curve, for the dashed 
curve Ub was fixed to zero, (b) Variation of the d-distribution at = 15 ± 5° with 
zi + Z2 and the result of the fit (solid line) yielding 5 mrad. The dashed line 

shows the fit at 825 mm for comparison. 
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Fig. 5. Sketch of the orbit distortion expected for a simultaneous excitation of the 
horizontal steerer SH5 and SH9. The beam position is measured by EDDA at the 
indicated position. 

with two fit parameters d and aa bearing the information on the beam. The 
total number of counts for a given ^^-interval A'^(0d) is used for the normal- 
ization and kept fixed in the fit. 

For a given (pd, the value of d is given by Eq. (4) and 



/ i?2 
(Td y C0S2 0dCT2^ + _^ 



(6) 



when the vertical size of the target is neglected (cf. Appendix A). Thus all 
16 spectra are fitted simultaneously by four fit parameters: Xj,, yb, , and 
(T^. For the implementation of the fitting procedure the minimizer MIGRAD 
(variable metric method) of the CERN package MINUIT [7] is used. Typical 
spectra are shown in Fig. 4 (a) for a 25 ms wide bin in time - corresponding to 
Ap ^ 29 MeV/c during acceleration of the COSY beam. Clearly, differences 
of 0.5 mm in can be resolved by the data. 



4.1.2 Test of Method 

As a test of the sensitivity of the method, data at fixed momentum of 3.4 GeV/c 
was acquired while the beam was moved horizontally across the target by ex- 
citing two horizontal steerer dipoles of the COSY-lattice. Between the steerer, 
separated by a phase advance of about tt, focusing quadrupoles are located, 
resulting in the geometry of the beam excursion as depicted in Fig. 5. The 
strength of the magnets (originally at 0%) was ramped within 100 ms to 
-1-80%, then, after 10 ms, they were ramped hnearly to -100% within 1000 ms 
and finally -after another 10 ms- brought back to 0% within 100 ms. 

The result of the fit (Fig. 6) shows the response of the beam: The beam position 
Xb follows closely the strength of the magnet. The vertical target position yj,, 
the beam width Sxb^FWHM, and the angular resolution 50 remain basically 
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Fig. 6. Result of the fit for flattop data at 3.4 GeV/c where around 3000-4000 ms 

a horizontal excursion of the beam was produced by exciting steerer magnets (cf. 
Fig. 5). The solid line shows the expected shape of the excursion with the amplitude 
and position in xi, adjusted to the data. 

unchanged. This proves the sensitivity of this method to the beam position. 

Another test is to take data taken at fixed energy and beam position and to 
observe the increase in beam size due to the emittance growth when the beam 
is "heated" by the target. The result of the fit is shown in Fig. 7. As expected 
only the beam size increases with time. All other parameters remain constant. 
It is interesting to see, that even the overshoot of the target - after is has been 
moved in vertically during the first 250 ms by a linear actuator and before 
it settles at the correct height - is resolved. As long as the target moves the 
vertex is then spread out in y and results in an increased angular width 50 in 
the fit. 



In all figures only statistical errors are given as they are calculated from the 
error-matrix of the nonlinear x^-fit. However, systematic errors arise mainly 
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Fig. 7. Result of the fit for flattop data at 2.25 GeV/c. The fiber target starts to 
move at ms and reaches the beam position around 250 ms. 

due to strong correlations of fit parameters. In Fig. 7 the anti-correlations of 
50 and Sxb^FWHM introduced by Eq. (6) is clearly visible. Since RS(/)/V2 
3.6 mm the angular resolution in (j) makes it difficult to determine beam sizes 
smaller than 3 mm accurately. Then, the cos^ (f)fi modulation of becomes 
very small and statistical fluctuations start to mislead the fit. 



4.2 Tilted Beam 



In the previous section it has been assumed that the beam is parallel to the 
symmetry axis of the detector. This is not necessarily true. In general, the 
beam may be tilted with respect to the detector both horizontally and ver- 
tically, described by the parameters and y^. Fig. 8 shows that a tilted, 
non-displaced beam mocks up a beam displacement in the analysis neglecting 
x'f, and y'l,: since only the azimuthal angles are considered, we reconstruct the 
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Fig. 8. Geometry for a tilted beam (axis z tilted by x^, y'^ as explained in the text. 

event in a plane parallel to the detector z-axis, defined by the two points of 
interception with the EDDA detector (dark shaded plane in Fig. 8). The reac- 
tion plane (light plane in Fig. 8) is defined by the same two points but must 
be coplanar to the beam. Consider a pp-elastic event where the two protons 
are emitted up-down: if the beam angle x'^ is positive the two protons tracks 
will be displaced to the left with respect to their origin. The displacement 
depends on the distance in z they travel until they hit the EDDA bar layer at 
radius R. As evinced in Fig. 8 the line through the two points of interception 
projected into the x-y plane looks as if the beam were displaced to positive x. 

In Appendix A it is derived that the measured value of d of a tilted beam is 
given by: 

d = {x^^ ^^-^x'f,) cos (pd + ivb + ^^-^y'b) sin (f)d (7) 

and reduces to Eq. (4) for = = 0. zi and Z2 are the 2;-coordinate of the 
point where the tracks of the proton intercept the EDDA bar layer. Evidently, 
the parameters Xb and x'j, (as i/b and yl) are strongly correlated and can only 
be distinguished by considering events with different zi + Z2- For pp-elastic 
scattering the correlation of the two z values as a function of beam kinetic 
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Fig. 9. Result of the fit for flattop data at 3.4 GeV/c. The solid symbols are the 

result of the fit when the beam angles .x^ and y'^ are included in the fit. The open 
symbols are the data as in Fig. 6 and the solid lines the expected shape of the 
excursion. 

energy Tp is given by Eq. (1) 

-1^2 = i?^e'.„.. = i?^(l + ^) (8) 

where rup is the proton mass. The values covered by zi + Z2 range from 
2-R7c.m. ~ 400 . . . 500 mm for symmetric events [zi — Z2) to 1100 mm for 
asymmetric events (^i = ^max = 1 m, see Fig. 1). 

4-2.1 Experimental Input 

For each ^^-bin the data was further subdivided into 14 classes with respect 
to their zi + Z2 values, ranging from 400 to 1100 mm, each 50 mm wide. For a 
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given energy only those classes were included in the fit where the kinematics 
allowed for the respective Zi + Z2 value. Usually 8-13 classes arc useful, and 
for each class 16 d distributions for the different 0^ ranges are sorted. The 
resulting 128-208 spectra are fitted simultaneously with six parameters x^, 1/6, 
(7^, (Tx^-i a^ft, and y^. An example of these spectra is shown in Fig. 4 (b) where 
the shift in d with zi-\-Z2, containing the information on beam angles is clearly 
visible through comparison to the result for zi-\- Z2^ 825 mm. 

The result of the fit (same data as in Fig. 6) is shown as the sohd symbols in 
Fig. 9. The change in beam angle and the beam position have the same sign 
as one expects for the geometry shown in Fig. 5. The total change in angle 
during the bump (4 mrad) and the range of the beam position shift (8 mm) 
yield a distance of the origin of the beam rotation at -2 m with respect to 
the target location which agrees rather well with the distance of the steerer 
magnet SH5 to the fiber target (1.75 m). The fact that the change in X}, and 
have the same sign indicates that the method is indeed sensitive to changes in 
beam angles, since Eq. (7) implies that fluctuations of these parameters due 
to statistics would be anti- correlated m. the fit. Again errors shown are purely 
statistical, and the observed fiuctuations around the expected excursions are 
now much larger, due to the strong anti-correlation of beam position and 
angle in the fit. These fluctuations can be only be damped by improving on 
the statistical precision of the data. Summarizing, possible beam tilts x'^ and 
y'^, should always be taken into account since the results for the beam positions 
and yi, depend on x'^ and y'^ (compare Figs. 6 and 9). 



5 Results and Discussion 

The method of reconstructing the first and second moment of the vertex distri- 
bution has been applied routinely for data analysis of the EDDA experiment. 
In Fig. 10 we show an example of a measurement taken during acceleration of 
the COSY beam and the flattop at 3.3 GeV/c. For comparison both the beam 
momentum and the response of the secondary electron monitor (SEM) as a 
measure of the luminosity are shown. In this example the beam position could 
be held rather stable throughout the ramp, some shifts in the beam angle a;^ is 
seen, especially in the transition to the flattop. Most strikingly, the rise in the 
SEM signal coincides with a dramatic increase in the horizontal beam width. 
Beam profile measurements showed a correlated decrease of the vertical beam 
width. Thus, the increase in luminosity on the horizontally mounted EDDA 
fiber target is due to a change of shape of the beam ellipse in the x-y plane 
at the target location, probably caused by a resonant coupling of the vertical 
and horizontal phase-space. The results show that including the beam angles 
in the fit is very important to reliably extract the beam position, besides the 
per degree of freedom is considerably reduced. The achieved resolution is 
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Fig. 10. Fit parameter as a function of time for data acquired during acceleration 
and the flattop. Also shown is the momentum of the beam and the signal of a 
luminosity monitor, the secondary electron monitor (SEM). Solid (open) symbols 
and solid (dashed) lines show the result of the fit when x|, and were (were not) 
included as fit parameters. 
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better than 0.5 mm for beam position and width and about 0.5 mrad for beam 
angles. 

Statistics obtained at typical running conditions of the EDDA experiment 
are sufficient to obtain this information with 25 ms resolution, corresponding 
to less than 30 MeV/c in momentum bins during acceleration with dp/dt ~ 
1.15 (GeV/c)/s ramping speed. Limitation of this method are, that displace- 
ments (<10 mm) and beam angles (<10 mrad) must be sufficiently small, as 
it was routinely achieved at COSY. Furthermore, beam widths below 2.5 mm 
FWHM cannot be determined very accurately, due to the limited azimuthal 
resolution. 



6 Summciry 

Exploiting the fact that pp-elastic events are coplanar, a measurement of 
the position where the particles intercept a cylindrical scintillator hodoscope, 
yields information on the position, width and angle of the impinging beam, 
without individual tracing of the particle tracks. The horizontal beam position 
and the vertical target position can be obtained with an accuracy of about 
1 mm, beam angles to about 1 mrad and the beam width to better than 1 mm 
FWHM (provided the beam is at least 2.5 mm wide). This information is ob- 
tained for time intervals of 25 ms which is sufficient to follow the smooth drift 
of a stored, accelerated proton beam in a synchrotron like COSY. 

We gratefully acknowledge excellent beam preparation by the COSY operating 
staff, including the setup of the beam excursions used for consistency tests. 



A Derivation of the Formalism 



A.l Coordinate Systems 

When the beam is both displaced and tilted with respect to the detector it 
is convenient to introduce a "beam coordinate system" BC {x,y,z) where the 
reaction vertex is at (0,0,0) and the beam parallel to the 5-axis. Transfor- 
mations from BC to the "detector coordinate system" DC {x, y, z) are then 
described by a rotation by an orthogonal matrix Ai followed by a translation 
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to the vertex position {xb,yb,0) (we assume the vertex is at z—0) 



X Xb 

y 

z 



Vb 



+ M 



X 



(A.l) 



The transformation matrix can be derived by expressing the unit vectors e^, 
By, and Sz of the BC system in the DC space for the special case Xb=yh=^'- 

• the unit vector in z is given by the beam parameter x'j, = dx/ds and = 
dy/ds. If one defines angles 9^ and 6y via tan 6*3; = x' and tan6'j^ = y' one 
obtains = (sin^j,, sin^,^, fc), where k is fixed by normalization 



k = \Jl- sin^ 6*2; - sin^ 9y 



(A.2) 



• for convenience we choose to be in the x-z plane. Orthogonality to and 
unit length fix to be {k/ cos9y, 0, — sin6'a;/ cos6'y). 

• ey is then given by the cross product x e^. 

The transformation matrix is thus given by: 



M 



k/ cos 6y 




■ sin 6x tan 6y sin 6^ 



^ — sin 9x/ cos 9y — 



cos 9y 

A; tan 



sin 9y 
k 



(A.3) 



The inverse transformation is obtained by inverting M.. Since A4 is orthogonal 
we simply use the transposed matrix: 



y 
V) 



I \ 

X — Xb 

y-Vb 

z 



(A.4) 



A.2 Dependence of d on Beam Parameters 



The quantity d (as explained in Section 1) is given by 



d — —Rcos 



>i — <P2 



(A.5) 
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as the average d for many events, so that smearing due to the finite detector 
resolution is averaged out. 

Lets now consider a perfect proton-proton scattering event with azimuthal 
angles 0j and polar angles 9i in the BC system. The unit vectors of their 
direction are given by 



Pi = 



sm (Ji cos (pi 
sin 9i sin 0j 
cos 9j 



1,2 



(A.6) 



We now calculate the point of interception of the two prongs with the cylin- 
drical bar layer of the ED DA detector at radius R in the DC system viz 



K cos (pi 
R sin (pi 

Zi 



Vb 



+ Xi Mpi 



; ^ = 1,2 



(A.7) 



by determining the three parameters Zi, (pi (describing the point where the 
bar layer is hit) and Aj. The result can not be obtained analytically, however, 
for some selected values for Xb, i/b, x'b, and y'b the general dependence can be 
extracted. Therefore, we first approximate the matrix for small angles x'b 
and y'b with the result 







Vb 



-y'b 



1 '2 '2 



^-Vb' 



Xi, 



y'b 



\ 



1 '2 '2 

1 - ^6 - yb 



(A.8) 



and consider the following cases: 
x'b ^y'b^O - 

Here, the beam axis and the detector symmetry axis are parallel, A4 is the 
unit matrix and we can write : 



R sin (pi — yb = Zi sin (pi tan 9i 
R cos (pi — Xb = Zi cos (pi tan 9i. 



(A.9) 



Dividing both equations and using coplanarity (tan^i = tan02) we can 
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write after some algebra exploiting trigonometric relations: 

01-02 01 + 02 , .01 + 02 / A 1 

i^cos ^Xftcos hyfesm (A. 10) 

which is the same as Eq. (4) when applying the sign convention for d and 
0d as explained in Section 4.1. 

Xb = yb = y'b = 0, 4 7^ : 
Here we derive: 

R cot 0i = ^Jl — x'^ cot j + 



Xu 



x'^ cot 0j + 



1 - % 



tan sin 0j . . s 

— , (A.ll) 



-Rsin0j tan sin 0j 

eliminating 6^ yields 

yr^cot0. - = (1 - 2x,') cot0.. (A.12) 

^ it sm 0j 

Coplanarity forces cot 0i = cot 02 and after some more algebra we obtain: 

01-02 f ^1 + ^2 01 + 02 

R cos = , < — - — cos 



2 2 2 



b ^ (A.13) 

^1-2:2 . 01 + 02 ^ 01 - 02 1 

— - — sm cot > 



Keeping only terms linear in x'^ and considering that ^i^^ ~ if for elastic 
scattering such that the cotangent is close to zero we may write: 

01-02 ,Zi + Z2 01 + 02 (K^A\ 

Rcos — Xh cos . A. 14 

2 2 2 ^ ^ 

An analog expression is straightforward to derive for x'j, = and 7^ 0. 

Taking the definition of d and 0^ from Section 4.1 we may now guess the 
dependence of d on Xb, y^, x^, and y^: 

if Z1+Z2 ,\ , / -21+-22 /\ . , /A 

d = ( Xb H — Xb 1 cos 0d + [Vb H ^ — Vb j sm 0d- (A.15) 



This result has been checked numerically to be accurate within 0.2 mm for 
reasonable values of beam parameters (|a;;,|, \yb\ < 10 mm,|a;^|, ||/^| < 10 mrad). 
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A. 3 Dependence of a a on Beam Parameters 



The width of the distribution in d at a given angle (f)^ is derived simply by 
using error propagation from Eq. (A. 15). 



dd V ( dd \^ ( dd \^ ( dd 



^'^^si,-''^} n^^N Ms^'-'J 

dd \ ( dd \ f dd ^ 



^[ar."-^) ^[a;;^") ^[w/*') '^'^'^ 



We derive 



2 

+ (xfe^ cos^ + sin2 0^) ^ (A. 17) 

+ I l/fe cos (f)d - Xb sin (f)d + ^ ^ ^ (t/b cos 0^ - x'^ sin 0^) | o-^^ , 

where the error in zi and Z2 was assumed to be the same. Here, the error of 
(pd is given by the detector resolution in the azimuth 

= ^ (A.18) 



Of all these terms most are negligible. Considering that x^, and yh are of the 
order mm, xj,, and of the order mrad, Zi of the order 400 mm, cr^ of the 
order 10 mm, and a^^ ft; 10 mrad, the last two terms are suppressed by two 
orders of magnitude with respect to the leading terms proportional to and 

<■ 

The width in Xh and x'^ are related by the amplitude-function oi the COSY- 
ring at the target location 

= ^ (A.19) 



provided we are at the beam waist. At the EDDA target location /3x and 
is of the order of some meters, such that we expect a^'^ and ay'^ to be less 
than a mrad. Therefore the terms proportional to crS and a^, will be a small 
correction to a'i^ and a"?.^ and we write 

'*'o yb 

~ 0-1 cos^ 0d + al sin^ 0^. (A.20) 
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This is only the dependence which enters through the right-hand-side of Eq. (A. 15). 
The quantity d itself is computed from the azimuthal angles which introduces 
another error through Eq. (3). This modifies Eq. (A. 20) to: 

a/ «i cos^ + aj^ sin^ + — ^ (A.21) 



where we have used that sin ^iz^a. ^ coplanar events. This expression, 
however, is not suitable for a fit. The three parameters ax^, Cy^, and cr^ are 
highly correlated, i.e. a combination of two can mock up the influence of the 
third. This is easy to see, if one multiplies the last term with 1 = sin^ 0^ + 
cos^ (pd- The assumption of a perfect angular resolution {a^ = 0) has the same 
effect as an increased size of the vertex area both in x and y. For a fiber target 
as it is in use with EDDA the vertex distribution is very narrow in y (some 
fim) so that without loss of accuracy may be set to zero and we obtain: 

a/ = <cos^0rf+^ (A.22) 

with two linearly independent parameters and a^. Alternatively, if the 
angular resolution is known from other sources, ax^ and may be extracted 
by this method. 
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